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Abstract 



We consider quantum theory of fields <j> defined on a D dimensional 
manifold (bulk) with an interaction V(<j>) concentrated on a d < D dimen- 
sional surface (brane). Such a quantum field theory can be less singular 
than the one in d dimensions with an interaction V ((/>). It is shown that 
, scaling properties of fields on the brane are different from the ones in the 

bulk. We discuss as an example fields on de Sitter space. 

1 Introduction 

Models with an interaction concentrated on a d < D dimensional submani- 
y\ ' fold (brane) of a D dimensional manifold (bulk) are interesting for high energy 

^ . physics as well as for statistical physics. In the first case we consider the visible 

universe as a submanifold ( a brane [l][2])of a higher dimensional space. Field 
theoretic models with an interaction on the boundary come also from string 
theory [3]. In ref. [3] the eleven dimensional gravity is interacting with ten di- 
mensional gauge fields living on the boundary. In statistical physics we may 
consider materials with a boundary and an interaction of some constituents 
placed on the boundary [4] [5] [6] [7] . It is an experimental fact [8] that correla- 
tion functions of field variables depending on the boundary points have critical 
exponents different from the bulk correlation functions. 

In this paper we discuss field theoretic models with an interaction on the 
brane. We concentrate on the scalar field but some methods and results can 
be generalized to gravitational and gauge field interactions. We begin with the 
free propagator. We admit any boundary condition preserving the symmetries 
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of the free Lagrangian. We show that a differential operator which is singular 
close to the brane has the Green function which is more regular on the brane 
than the one for operators with constant coefficients. Subsequently, we discuss 
models with an interaction concentrated on the brane. The functional measure 
is defined [9] by its covariance (the Green function) and its mean. The mean 
breaks symmetries of the classical action. We average over mean values in order 
to preserve the symmetries. As a consequence of the more regular behaviour of 
the Green function the model with an interaction concentrated on the brane has 
milder ultraviolet divergencies. We give examples of nonrenormalizable theories 
in the bulk which become superrenormalizablc when restricted to the boundary. 
We work mainly with the imaginary time version of quantum field theory. In 
the last section we discuss the scattering theory. The free particle is treated 
as a packet of waves on a curved manifold. We calculate the scattering matrix 
of such particles resulting from the V{<j>) interaction concentrated either on the 
boundary or at the time z = (a kick at a fixed moment). 



2 Green functions on a boundary 

We consider a D = d + m dimensional Ricmannian manifold of the warped 
form M g — M m x g R d [10] with the boundary R d whose metric close to the 
boundary takes the form 

ds 2 = G AB (X)dX A dX B = g lLV { i) )dx»dx v + g jk {y)dy j dy k (1) 

where X — (y,x) are local coordinates on M g , gjk is the Riemannian metric 
induced on M m and g^ : M m — » R d is a positive definite d x d matrix function 
defined on M. m . The action for the free field <j> reads 

Wo = J dxVGG AB dA(t>d B <l> (2) 

The free (Euclidean) quantum field can be defined as the one whose propagator 
is determined by the Green function 

-AG = d A G AB VGd B G = S (3) 

where G = det Gab- In the metric (1) eq.(3) can be expressed as 

A + ^ k (y)VG(y)^)g = 8 (4) 

The solution of eq.(3) is not unique. If Q' is another solution of eq.(3) then 
Q' = Q + 1Z where 1Z is a solution of the equation 

ATI = (5) 
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We can determine Q unambiguously imposing some additional requirements, 
e.g., requiring that Q = on the boundary or that Q be scale invariant on the 
boundary. 

We assume that the metric is scale invariant 

g^(\y) = ^ 2a 9^(y) (6) 

and 

g ok {\y) = \ 2f3 g ]k (y) (7) 

It follows that 

G(Xy) = \ 2da+2mf> G(y) 
Let us consider G\ tP (x, y. x', y') = G(px, Ay; px' , Ay'). It satisfies the equation 



X m pd (x-^+™P+^g^{y)VG{y) ^ ^ 

+X mp +da -2,_a_ ig3k{y)V Q {y) _a_^ g = §{x x , )s{y yl) 



(8) 



Eq.(8) is identical with eq.(4) but expressed in rescaled coordinates (the scale 
invariance (6)- (7) has been applied). Let us choose p such that the scale factors 
in the two terms in eq.(8) are equal 

x -2u+mP+da ^-2 _ ynP+du-2f}-2 ^g-j 

Hence, 

p = x i- a +p (1Q) 

Then, it follows that 

Or MVGfo)^ + ^ k (y)VG iy) ^)g = 6 (11) 

where 

G(y, x- y', x') = A^^+^Ay, A 1 "^*; Ay', A 1 "^*' ) (12) 

Q and Q satisfy the same equation (4). Hence, for the scale invariant solution Q = 
Q (there are many solutions of eq.(4) but the scale invariant solution is unique). 
We are interested in the boundary Green function (5(0, x; 0, x') = Qe{\x — x'\). 
Qe depends solely on the Euclidean distance |x — x'\ as a consequence of the 
translational invariance of eq.(4) and its rotational invariance when y = y' = 0. 
Choosing in eq.(12) 

A = \x-x'\~^+e (13) 

we obtain 

g E (\x-x'\) = K\x-x'\-( D - 2 ^ (14) 

where 

a = T^Tp (15) 
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and K(a,f3,D) = Q(0, e; 0.0) (here e € R d is an arbitrary vector such that 
|e| = 1). Some authors [?][?] [11] in de Sitter case require that Q = at y = 0. In 
such a case K = in eq.(14) and our result is trivial. However, in other models 
of QFT on dc Sitter space [12] [13] [14] [15] [16] Q docs not vanish at y = 0. 
It will be useful to rewrite eq.(4) in the momentum space. Let 

9(y, x; y', x') = (2^)-" 1 j dpcMw(x - x'))G(p; y, y') (16) 

Then, Q satisfies the equation 

( - P,Pug^{y)VG(y) + A^ (j/) VG( y )^-)e = 6(y y') (17) 

Let Gx,-y(p;y,y') = G{lP; Ay, Ay'). Then, is a solution of the equation 

+x -23 +mf s +da _^ gjk{y)V G {y) _a_ w ^ = 6{y _ yl) 

We choose 

Then, we obtain (similarly as in eq.(12)) from the uniqueness of the solution of 
eq.(17) 

G(p; y, y') = X^-^+^+^giX-^+y- \y, Ay') (19) 

Setting 

A = \p\T=i+i3 

and y = y' = in eq.(19) leads to the result 

G(p;0,0) = K\p\- 2 " (20) 

where 

-2. = + /?)+*» (21) 

l + /3-a y ' 

and K(a, (i, D) = G(e: 0, 0) ( e is an arbitrary unit vector in R d ). 

As an example of an application of eqs.(20)-(21) we could consider de Sitter 
space with the metric (it is of the type (l);this is Euclidean AdSc [18], the 
notion of the boundary in AdS and its relation to Euclidean AdS is discussed 
in [17] [18]) 

ds 2 = dt 2 + cx_p(2Ht)(dxl + .... + dx 2 D _ 1 ) = y~ 2 (dy 2 + dx\ + .... + dx 2 D _ 1 ) (22) 
where y = H^ 1 cxp(-Ht) . 
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Then, from eqs.(20)-(21) 

G(p;0,0) = K\p\- D+1 (23) 

The Fourier transform of Q in eq.(23) is infrared divergent. It can be defined as 
a distribution on a set of functions vanishing at p = 0. Then, 

g E {\x-x'\) = -K\n\x-x'\ (24) 

Eq.(24) gives the form of the Green function for a = in eq.(14). 

Generalizing eq.(4) we could consider a system of equations for a tensorial 
Green function Grn 

(25) 

here a sum over repeated indices £ is assumed, the Sra function is the product 
of the usual 5 function of eq.(8) and a scale invariant tensor with indices Tfl. 
We assume that the coefficients have the following scaling properties 

a"" (Ay) = \ {d - 2)a+mf3 a^(y) (26) 

V k (Xy) = \ ia+ ^ n - 2 ^V k {y) (27) 

A^(\y) = \( m ~ 2 ^ +da ~ 1 A(y) (28) 

B° (Ay) = \(™-i)0+(d-i)a-i B j( y j ( 29 ) 

and 

V(Xy) = \{™-W+da-2 V ^ 

Then, repeating our scaling arguments of this section we could derive the results 
(14)-(15) for Gm and (20)-(21) for its Fourier transform. Such results may be 
applicable to propagators describing an interaction of gravity on the bulk and 
on the brane with gauge fields on the brane as in [3][1]. 



3 DeWitt expansion 

Let us discuss now scaling properties of the Green functions from the point of 
view of DeWitt expansion [19] and the Hadamard representation of the Green 
functions [20]. According to the DeWitt suggestion we can solve the equation 
(for a non-negative operator .4) 

-AG = 5 (30) 

by means of the heat kernel which is the fundamental solution of the heat 
equation 

—^-JC T = AJC T (31) 
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Namely, 

9=1 drlC T (32) 
Jo 

Then, DeWitt makes the assumption 

Kr(X, X') = (2ttt)-* exp(-^%^)A(r, X, X 1 ) (33) 

zr 

where a(X,X') is the square of the geodesic distance between X and X' . A 
has a Taylor expansion in powers of r for a manifold without boundary (powers 
of i/t may appear if the manifold has a boundary) . Performing the integral 
over t in eq.(32) we obtain the short distance expansion of Q(X,X'). 

We do not know the formula for a in general. We can explain the expansion 
(33) in the hyperbolic case (22). Then, 

ch(a(X, X')) = 1 + {2yy')-\{x x'f + (y - y'f) (34) 

In the hyperbolic space the heat kernel is known exactly [21]. The DeWitt 
expansion holds true (A has an expansion in powers of er; no r dependence). 
We can see that 

(t(A,A') _1 ^0 (35) 

when y — > (t — > oo in cq.(22))and x ^ x' . It follows from eq.(32) that if 
D = 2k + 1 is odd then 

g E (x-x') = (36) 

and if D = 2k is even then 

Ge{x-x') = -Kln\x-x'\ (37) 
The conclusion (37) coincides with eq.(24). 



4 Interacting fields on the boundary 

In the model (2) we consider an interaction V concentrated on the boundary 
y = (we do not treat here a coordinate independent geometric description of 
the boundary but restrict ourselves to the model (1)) 

W = W + Wj = J dxVGG AB d A <j}d B (t> + J dX\/G5(y)V(<f>) (38) 

We define the functional measure 

dfi(<j>) = Z^V^cxpi-W) = Z^dnoi^expi-W!) (39) 
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where the Gaussian measure is 



diMj{(j)) = D0exp(-W o ) 

The partition function 

Z = j d/i exp(-W» 

determines a normalization factor. The Gaussian measure is defined [9] by the 
mean 

J d^)4>(x) = (cj>(X)) 

and the covariance 

g(x,x') = J d^mx) (<t>(x)))Mx') - mx'))) (40) 

In the papers on AdS-CFT correspondence [18] [24] [17] [?][?] the choice is made 
Q(X,X') = Gd(X, X 1 ) where Go is the Dirichlet Green function (vanishing on 
the boundary) and (<j>(x)) — ip(x) where ip{ x ) is the solution of the equation 

Atp = 

with a fixed boundary condition. However, the choice of the boundary field 
ip 7^ breaks the rotational and translational invariance in the x variables 
present in the classical action (38)with the metric (1). The approach with the 
classical boundary field and the Dirichlet boundary condition leads to a different 
quantum field theory than the one developed in refs.[12][16][13][14][22] (see also 
[15] [23]). In our approach the QFT is determined by the choice of the Green 
function Q (40) ( we set the mean (</>) = 0). We choose the Green function Q 
which has the symmetries of the action Wo (2). Hence, the functional measure 
(39) will have the symmetries of the action (38). 

We can show that our approach is equivalent to a quantization with a given 
classical solution ip if subsequently an average over all such solutions is per- 
formed. We assume that in the sense of bilinear forms 

G>G D (41) 

Such an inequality follows from the maximum principle for elliptic operators 
[25] [26]. The inequality (41) holds true also in the non-elliptic cases discussed 
in our earlier papers [27] [28]. Using eq.(41) we may write 

G{x,x') = G D (x,x') + g B {x,x') 

where Qb is a non-negative bilinear form and Qb (0, x; 0, x') =Qe[x — x'). Then, 
[9] [29] 

J duoi^cxpi-WT^Ficf)) = J dii D (<t> D )diiB(<f>B)exp(-Wi(<f> D +<f>B))F(4> D +4> B ) 
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where 4>d is a random field with the covariance Gd and 4>b is the random field 
with the covariance Qb- Let us note that because Gd as well as Q satisfy the 
same equation (2) then their difference satisfies the equation 

A(x)g B (x, x') = A(x')g B (x, x 1 ) = o (42) 

We can solve eq.(42) with the given boundary condition Q E 

g B (X,X') = J dx b ^g J dx b ^gV(X,x b )V(X',x b )G E (x b -x' b ) (43) 

where V is the Green function solving the Dirichlet problem (the boundary to 
bulk propagator). It follows that 

<j> B {X) = J dx b JgV{X,x b )3>{x b ) (44) 

where $ is the Gaussian random field defined on the boundary with the covari- 



ance 



($(x)$(a^)) - Q E {x - x') (45) 

In the case of the de Sitter space the solution of the Dirichlet boundary problem 
(44) can be expressed in the form 

4>{X)=yi dpexp(ipx)\p\%K 4 (\p\y)$(p) (46) 



where $ is the Fourier transform of $ and K v is the modified Bessel function 
of order v. From eq.(23)[27] 

(<S>(pW(p'))=5(p-p>)\p\- d 

The Schwinger functions are defined as moments of the measure \x 

^(x 1 )...4ix k )) = z- 1 f d/i(M(Xi) 4>(Xk) 



We calculate these Schwinger functions in the N-th order of the perturbation 
expansion 

{4>{x 1 )....<t>{x k )) N = 

Z" 1 / dX{....dX' N 6(y[)....S(y' N ) (47) 
Jd^oHXi) 4>(X k )V(4,(X{)) V{<j>{X' N )) 

If V(4>) is a normal-ordered polynomial of order r then the Schwinger functions 
of order N are expressed by a product of at most k Green functions Q (y, x; 0, x') 
and at most (rN)\ Green functions Q(0,Xj-,0,x' k ). It follows that if the Green 
functions at y = y' = are sufficiently regular (depending on a and (3 in 
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eqs.(14)-(15)) and the integration over x in the interaction Wi is restricted to 
a finite volume A then the Schwinger functions (47) for non-coinciding points 
Xj = (0,Xj) are finite and non-zero . In fact, the amputated Schwinger func- 
tions (when the propagators corresponding to the lines connecting the external 
points are removed) coincide with the ones for the V((f>) theory in R d calculated 
with the propagator (14). If in eq.(47)X, = (yj,Xj) and we set all yj — then 
the resulting perturbative quantum field theory coincides with the V(<f>) theory 
where the conventional propagators are replaced with Q E of eq.(14). In partic- 
ular, the partition function Z in eq.(40) can be finite and non-zero even though 
the bulk theory is non-renormalizablc. Note, that this result is a consequence 
of the singularity of the metric on the boundary. Adding the boundary inter- 
action Wi to the regular free part Wo would lead to the theory with the same 
singularity as the free field theory in d + 1 dimensions. 
Let us write the propagator (14) in the form 

Q E = K\x-x'\- d+p (48) 

then the effective field theory on the boundary y = could be represented by 
an equivalent functional measure fi of the form (39) with 



W 



= W + W I = c J dx(t>(-A)2<t> + J dxV{(j>) (49) 



where A is the Euclidean d-dimensional Laplacian. If the theory (49) is to be 
scale invariant (up to the log-terms coming from the renormalization) 

4>{x) ~ \ u 4>{\x) (50) 

then from the propagator (48) it follows that 

d-p = 2u (51) 

Hence, the interaction 

' dxV((f>) = n I dx(f> r (52) 



is scale invariant if the order r of the interaction is related to v (and by eq.(51) 
to p) 

v=- (53) 
r 

From eqs.(51) and (53) it follows that 

p = d(l--) (54) 
r 

If p = d then from eq.(54) we obtain r = oo. This case corresponds to the 
exponential potential V(cj>) = Kexp(^). In general, we obtain simple fractional 
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scaling dimensions (54) determined uniquely by the natural numbers d and r. 
From eq.(15) and eq.(54) 

2,= (d-l )T l±A_ ( 55 ) 

Hence, the geometry of the surface imposes conditions on the scaling exponents 
and the form of the scale invariant interaction. From eq.(53) r = ^ hence 
the order of the scale invariant interaction is also determined by the surface 
geometry 

r - (d-l)(l + /9) (56) 

If (3 = — 1 then f = and p = d. In such a case from the scale invariant fields 
(24) (of dimension zero) we can form conformal fields of higher dimensions by 
means of exponential functions in a similar way as in two dimensional conformal 
field theory [31] . 

Let us note that the theory with an interaction on the boundary can be 
considered as a scaling limit of the one with an interaction in the bulk. For this 
purpose define V\((j>(y, x)) = V((j)(Xy, x)), calculate the Schwinger functions (47) 
perturbatively and at the end take the limit A — > 0. 

5 M d+1 — > R d reduction 

If m = 1 then we may change coordinates in eq.(4) introducing a new coordinate 
z instead of y in such a way that 

f z =9 DD {y)VG{y) (57) 



Then, eq.(4) reads 



[g^{y)G{y)g^{y) — — +di)g = g DD VGS(y-y')S(x-x') = 6(z-z')6(x-x') 

(58) 

The action (38) takes the form 



W = J dydxVG(g DD d v <pd y <P + g^d^d v <P + Ko8(y)V(4>j) (5Q) 
= / dzdxd z 4>d z 4> + J dzdxg DD Gg>* u d^d„(j) + K J dzdxVGS{z)V \4>) 

In scale invariant models G(0) is either zero or infinite. Then, we must renor- 
malize the interaction defining kq = koVG(0). We have studied the behaviour 
of the Green functions Q of eq.(58) in refs.[27] [28]. Assuming that 

g^g DD G~ 5^\z\ 2 ~< (60) 
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for a small z we have shown (this is the same result as the one in eq.(14)) that 
for small \x — x'\ 

g E (x-x')~\x-x'\- d+ ^ (61) 

As discussed in sec. 4 the theory at z = is the same as the one with the 
interaction V and the propagator (48). Hence, depending on the value of 7 > 
— 1 + ^ the boundary field theory can be much more regular than the bulk field 
theory . When 7 = ~1 + 3 ( as in the case of de Sitter space) then the model 
(59) becomes superrenormalizable for polynomial interactions. 

It is interesting to consider d = 2 and p — 2 with the exponential interaction 
exp <f> on the boundary. Then, the boundary correlation functions are conformal 
invariant (the Liouville model [30]) whereas the bulk correlation functions are 
not renormalizable. 

The other interesting case is the boundary of the Euclidean (anti)de Sitter 
space [18](eqs.(22)-(24))with an exponential interaction. Again the boundary 
correlation functions of exponentials can be conformal invariant [31] whereas 
the bulk correlation functions are non-renormalizable. 

6 Semi-infinite statistical systems 

In this section we relate the model of sec. 4 to some models of statistical physics 
[4] [5] [6] . Let us consider Ginzburg-Landau Lagrangians with a scalar field <j){X) 
where X e R D , X = (z,x) with x G R d and z > 0. The boundary is at z = 0. 
We consider the action 

W{(f>) = JdzJ dx(c x Y%=i d^d^ + c 6(z) Y%=i d^d^cj) + c 3 d z (j)d z (j) 
+k 1 V 1 (4>) + koJ(*)W)) 

(62) 

The free propagator for a perturbation expansion around the Gaussian theory 
is determined by the equation 

d d 
(ci ]T fyfy + c S(z) 9^ + Ci d z d z ) Q = S (63) 

^L — 1 11— 1 

This model is related to the models of sees. 4 and 5. We have studied the 
propagator (63) in rcfs.[27][28]. If in the equation (59) = g^g DD G = 
t nen the scaling properties of a^ v are the same as the ones of S(z) 
(i.e., d(Xz) = \~ 1 5(z)) . In such a case the propagator (63) has the same short 
distance behaviour as the one for a^ v {y) — S^^l" 1 and b = 1 in eqs.(26)-(27) ( 
if Ci = then these propagators coincide). Now, in eq.(26) (d — 2)a + m/3 = — 1 
and (to — 2) [3 + da = 0. Hence, in eq.(14) <r(d — 1) = d — 2 and if c\ — then 
we have exactly 
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g E (\x-x'\) = K\x-x'\- d+2 (64) 

As a consequence the propagator in (d + 1) dimensions behaves as the one in 
d dimensions. In particular, if d + 1 = 3 then with the gradient term on the 
boundary (co ^ )we obtain the logarithmic short distance behaviour of bound- 
ary correlation functions. In such a case the model on the boundary (ki = 0) 
has the properties of the two-dimensional Euclidean field theory. If Vq ((f)) — <fi 6 
then we obtain a renormalizable field theory in the bulk which on the boundary 
reduces to the well-studied model of a tricritical phase transition . If d + 1 = 4 
and Vq(4>) — </> 4 then we have a typical superrenormalizable Landau- Ginzburg 
model. The model is superrenormalizable because the truncated Green func- 
tions are the same as in the threedimensional theory. The model in d + 1 = 5 
with the 4 interaction on the boundary is still renormalizable. In terms of 
the lattice approximation our semi-infinite model with c\ = describes spins 
whose interaction in the bulk is restricted to the lines perpendicular to the 
boundary. It is still surprising that although the system is D dimensional the 
correlation functions behave like the ones of an <i-dimcnsional system in spite 
of the opportunity for the interaction to spread into the D-th dimension. 



7 Scattering theory with an interaction concen- 
trated on the boundary 

We could approach the models of sees. 4-6 in the conventional Hamiltonian 
framework. We consider the scattering theory for the field equations 

A<f> = -5(z)V'(<f>) (65) 

We could treat z either as a spatial coordinate or as a time. Let us concentrate 
here on the latter interpretation. Then, for z — > — oo the interaction is switched 
off and <f> — > (p m where 

A4> m = (66) 

We consider here the real time z and the wave operator A = —d 2 — A 2 where A 2 
is a positive operator. We quantize the field <f> m and construct the interaction 
Hamiltonian Hj(z) = 5(z) J dxV(<f> m ). We can derive the S-matrix describing 
an interaction V(<j>) of the field <f> m by means of the conventional reduction 
formulas [16] [33]. We assume that the geometry is fixed and therefore the 
asymptotic states are constructed on a given gravitational background. As 
usual the S-matrix is determined by the time-ordered correlation functions 

(T(^{X 1 )...4{X n ))) = (T(ct> in (X 1 )...<p in (X n )eM-i j dzdx8{z)^gV{^ n )))) 

(67) 
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The expectation value on the rhs is calculated in the vacuum for the free field 
(f) m defined on the curved manifold (there may be many vacuum states and it 
belongs to the theory to choose one of them). In perturbation theory (in the 
Euclidean region)the correlation functions are expressed by the Green functions 
(47). The S-matrix is determined by r- functions 

S = 1 + E / dX 1 ....dX nT (X 1 ....X n ) : r n (X!) 4> ln (X n ) : +... 

n 

where 

t{X u X n ) = A(X 1 )....A(X n ){T(<j>(X 1 )...<p(X n y)) 

Hence, from eq.(47) the A operators cancel the external propagators and the r 
functions depend on Green functions defined at z = 0. 

Let us consider some examples. First, the case when in eq.(58) 

g DD Gg^(y) =S^\z\~ 1 (68) 

(this is an analogue of eq.(63)). Then 

A = \z\- 1 A-d 2 z (69) 

The solutions of eq.(69) have the form 

4? n (z,x) = ^zj d P {a{ P )H\ l \c\p\V~z)+a + {p)H^\C\p\V~z)) expiipx) (70) 

where is the Hankel function of order v [32] and C is a positive constant. 

A quantization of a and a + (as creation and annihilation operators) leads 
to a quantum scattering theory described by an expansion of the S-matrix in 
asymptotic fields <p m . Note, that at z — the Feynman propagator 

Q E (x -x') = K j dpcxp(ip{x - x'))\p\- 2 (71) 

for A in eq.(69) is equal to the one for a massless Euclidean free field in d 
dimensions. 

As a second example we consider a V interaction of particles in the hyperbolic 
space (22) (whose analytic continuation y — > iz gives de Sitter space). In such 
a case eq.(66) has the solution (at the real time y ~ z?) 

<f> in (y,x) =r J dp(a(p)H^(C\p\y) + a+(p)H^ (C\p\yj) exp(^) (72) 
where ui = ^. 

If the interaction is of the form 8{y)V then the propagators are logarithmic 
and the model will be superrenormalizable for (normal ordered) polynomial V. 
There will be no ultraviolet divergencies in the S-matrix. 
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8 Summary 



We find it as a remarkable property of some singular second order differential 
operators in D dimensions that the Green function restricted to the boundary 
can be less singular than \x — y\~ D+2 . When we put an interaction on the 
boundary then the regularity of the model can be extended to models with a 
non-trivial scattering. Models with an interaction on the boundary are certainly 
of physical relevance in statistical physics. It is still unclear whether such models 
are acceptable in high energy physics as , e.g., the brane theories. We suggest 
however that restricting an interaction to a brane can give a promising way of 
avoiding ultraviolet divergencies in quantum field theory. Our construction of 
the quantum field theory with a boundary ( which preserves the symmetries of 
the boundary) is different than the AdS-CFT approach. However, there appear 
some interesting relations between both approaches if we treat the boundary 
value (f> a as random (quantum) and average over <p Q . 
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